Quasiprobabilistic Interpretation of Weak measurements in Mesoscopic Junctions 
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The impossibility of measuring noncommuting quantum mechanical observables is one of the most 
fascinating consequences of the quantum mechanical postulates. Hence, to date the investigation of 
quantum measurement and projection is a fundamentally interesting topic. We propose to test the 
concept of weak measurement of noncommuting observables in mesoscopic transport experiments, 
using a quasiprobabilistic description. We derive an inequality for current correlators, which is 
satisfied by every classical probability but violated by high-frequency fourth-order cumulants in the 
quantum regime for experimentally feasible parameters. 



Every measurement in quantum mechanics is in prin- 
ciple described by the projection postulate However, 
in practice perfect projective detectors often do not ex- 
ist and the measurement encounters a finite error. This 
can be resolved by replacing the projection by Kraus 
operators defining a positive operator-valued measure 
(POVM)^^]. The Kraus operator can be continuously 
changed from projection - strong measurement (exact) 
- to almost identity operator - weak measurement (with 
huge random error). Effectively, a POVM means that we 
take detector's degrees of freedom as part of the consid- 
ered Hilbert space and make a projective measurement 
on the detector. Obviously, in that case the detector- 
system coupling defines the strength of the measurement 
of the system. The equivalence of a POVM and the pro- 
jective measurement follows from Naimark theorem 
The actual modeling of a detection scheme by POVM is 
a long-standing problem 0-0] • 

The interpretation of the results of a weak measure- 
ment can lead to paradoxes. For instance, if a weak 
measurement of A performed on the state p is followed 
by a projection B then the weak value can be defined 
b{A)p = Tr(Sip)/Tr(S/5) H. The unusual feature of 
the weak value is that it can exceed the spectrum of A, 
which obviously contradicts our classical intuition. The 
strange properties of weak measurements have been con- 
firmed experimentally in quantum optics j^, while ex- 
periments in solid state physics are proposed [lo| . The 
interpretation of recent experiments on current fluctua- 
tions in mesoscopic junctions in the quantum regime 11 1 
is impossible in terms of a usual probability [l2j . Instead, 
we proposed to consider a weak current measurement, 
which implies a large background noise, but avoids the 
paradox of a certain average square to become negative. 
The necessity of a weak measurement lead to corrections 
in the observed finite- frequency noise, however, a direct 
experimental proof of this scheme was not feasible. 

In this Letter, we first construct a general formula for 
a quantum quasiprobability, which does not depend on 
the details of the measurement apparatus and confirms 
the previously used formulas 12|. Second, we propose a 



scheme to test experimentally its negativity in frequency 
domain. To this end, we will derive a classical inequality 
for high-frequency current correlators of the form 

CL, < (a„ + 2ttSI/5uj) {C^,^, + 2^Sl,l5u') . (1) 

Here C^^, = {{I{uj)I{~uj)I{lo')I{-~uj'))) /to is a fourth- 
order correlator and — {{I{uj)I{—uj))) /to the 
frequency-dependent current noise, where to is the total 
(long) measurement time and Su) is the bandwidth of the 
detector {S and C are independent of to). Inequality ([T]) 
is satisfied by every classical stochastic process, but can 
be violated by high-frequency correlators in the quantum 
regime of a mesoscopic junction for experimentally acces- 
sible parameters. We believe our proposed violation of 
the classical inequality ([T]) can be realized with the ex- 



isting techniques [11|. This violation will be a proof of 



negative values of the quasiprobability. Although it is not 
necessary to explain the strange features of weak values 
[S], it offers an alternative test of nonclassicality similar 
to the Wigner function [l3|. Moreover, the quasiprob- 
abilistic interpretation can be easily generalized to an 
arbitrary sequence of measurements. This interpretation 
facilitates the transfer to mesoscopic junctions and we 
present an example, how the negativity of the quasiprob- 
ability can be proven in a tunnel contact and discuss the 
experimental feasibility. 

We will construct the quasiprobability by a deconvo- 
lution from a suitable POVM. The real parts of weak 
values can then be expressed as averages with respect to 
the quasiprobability. Let us begin with the basic prop- 
erties of a POVM. The Kraus operators K{A) for an 
observable described by A with continuous outcome A 
need only to satisfy / dAk'<{A)K{A) = i. The act of 
measurement on the state defined by the density matrix 
p results in the new state p{A) = K{A)pK^A). The new 
state yields a normalized and positive definite probability 
density p{A) = Trp(A). The procedure can be repeated 
recursively for an arbitrary sequence of (not necessarily 
commuting) operators Ai, . . . , An 14!], 



,An)=K{An)piAu...,Ar,-i)K\An) . (2) 
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The corresponding probabiUty density is given by 
p{Ai,...,An) = Tr p{Ai, . . . ,An). We now de- 
fine a family of Kraus operators, namely Kx{A) = 
(2A/7r)i/4exp(-A(i - A)'^). It is clear that A ^ oo 
should correspond to exact, strong, projective measure- 
ment, while A — > is a weak measurement and gives a 
large error. We also see that strong projection changes 
the state (by collapse), while A — >■ gives p{A) ~ p, 
and hence this case corresponds to the weak measure- 
ment. However, the repetition of the same measurement 
k times effectively means one measurement with A — > fcA 
so, with k ^ oo, even a weak coupling A <C 1 results 
in a strong measurement. For an arbitrary sequence of 
measurements, we can write the final density matrix as 
the convolution 



px{A) = J DA' gxiA')]lgkiAk - A'^) (3) 



with gk{A) 



\/2Xk/n. Here A = (Ai, . . . , A„) 



A — {Ai, . . . , An), and DA = dAi . . . dAn- The quasi- 
density matrix g is given recursively by 



SxiA) = I -e 



(4) 



with the initial density matrix g = p for n — 0. We 
can interpret g in ([3]) as some internal noise of the detec- 
tors which, in the limit A — > 0, should not influence the 
system. We define the quasiprobability gx — Tr gx and 
abbreviate g = go. In this limit (j4]) reduces to 

dx 



<>(A) ^ I ^- 



g{Ai 



, An- 



1 e 



(5) 

Note that go...o,A — g, so the last measurement does 
not need to be weak (it can be even a projection), and 
marginal distributions are consistent with absence of a 
measurement, / dAk g{A) = g{. . . , Ak-i, Ak+i, ■ . ■). In 
the case of commuting operators the quasiprobability re- 
duces to the usual probability g = poc- For Ai = x 
and A2 — "p with [x,p\ — ih we obtain the Wigner func- 
tion g{x,p) — g{p,x) — W{x,p) The definition 
preserves locality - for (A, pa) and {B,ph) acting in two 
separate Hilbert spaces we have p — paPb — 5> g{A, B) = 
ga[A)gh{B). The averages with respect to g are eas- 
ily calculated by means of the generating function ([5|). 
e.g. {A), = Trip, {AB), = Tr{i,B}p/2, {ABC), = 
TyC{B,{A,p}}/A for A ^ {A,B,CY This ordering of 
operators is called time symmetric [H, [3| • To relate the 
quasiprobability to weak values, we have to consider two 
measurements: A and B. The real part of the weak value 
is just the average Re b{A)p = {A)g\B with respect to 
the conditional quasiprobability g(A|i3) = g{A,B)/g{B). 
The complex weak values require a different interpreta- 
tion , which can also be generalized to sequential mea- 
surement 17 1 . 



We shall apply the above scheme to the measurement 
of current I{t) through a mesoscopic junction in a sta- 
tionary state. For a moment, we forget about quantum 
mechanics and recall basic properties of stochastic pro- 
cesses [3], applying them to SI ~ I — (/). It is con- 
venient to define the noise (second cumulant), S(a,/3) 
— 2nS{a + (3)Sa = {5 1 [a) 5 1 {(3)) , and the fourth cumu- 
lant C{a, /3, 7, 77) = 2'K5{a + (3 + ^ + ri)C{a, /3, 7,77) with 

C{a, /3, 7, v) - {5I{a)6I{(3)I{-f)6I{7j)) (6) 
-S{a, /3)^(7, v) ~ 5(a, j)S{P, 77) - 5(«, 77)5(7, /3). 

Here and throughout the text we use Latin arguments 
in time domain and Greek ones in frequency domain, 
related by a{uj) = J dt e*"*a(t). Note, that the delta 
function of the frequency sum has a cut-off of the order of 
the measuring time to (larger than all relevant timescales 
of the system) , which in some following expressions is a 
simple prefactor and docs not enter final conclusions. 
Let us define the fluctuating noise spectral density 
= J^^ SI{a)SI{~a)da with uj± = cjzLSuj/2, for which 
we obtain the average fluctuations 



{{SX^)^)/to= Co^^dadp + 27r Sida, 
{SX^SX^,)/to = da Ca^dp (7) 



where 5X = X {X) and Co,p = C{a, -a, /3, -l3). The 
intervals [a;_,w_|_] and are nonoverlapping. Con- 

sidering classical correlators of 6X at different frequencies 
we obtain the Cauchy-Bunyakovsky-Schwarz inequality 



B 



{SXJX^ 



{{SX^)^){{SX^ 



< 1. 



(8) 



If we choose e.g. < < cj^ < the correlators 
correspond to a low and high frequency measurement. 
Furthermore, assuming that S and C are constant within 
the bandwidth Slu, Slu', the inequality ([8]) takes the form 
([T]) mentioned in the introduction. It is interesting to 
note, that for frequency-independent (classical) noise the 
inequality is always satisfied. 

Turning to the quantum case we stress that continuous 
measurement cannot be strong as we would end up with 
the quantum Zeno effect and suppress the dynamics of 
the system completely [l9| . This follows from Eq. ^ 
for Heisenberg operators Ak ^ A{tk = fcA) and a finite 
tn — nA. In the limit ^ A^ — ?> 00 and n — >■ 00, g be- 
comes diagonal in the eigenbasis of A and freezes. Non- 
classical behavior of quantum correlations in the limit of 
weak measurement can also be shown using the Leggett- 
Garg inequality [13], which involves time- resolved sec- 
ond order correlations assuming that the observables are 
bounded. Our inequality is more general as we do not 
require bounded observables but we need fourth order 
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correlations in the frequency domain instead, which is 
more suitable for electric current measurements. 

We denote the measured current operator in the 
Heisenberg picture by I(t). Strong, projective measure- 
ment can be performed only if we are interested in the 
long-time limit. The finer the time resolution we want 
the weaker the measurement must be as [/(t), I{t')] 0. 
We define a generating functional of the quasiprobability 
in the weak measurement limit by 

g[I] = J D;^e-*^>^(*)^(*)^* (9) 
X Tr \Te-l' «x(t)/(t)dt/2 pf^J txit)i{t)dt/2 

with T (T) denoting (anti)time ordering. This represents 
a straightforward generalization of the generating func- 
tion obtained from the trace of ([S]), in which the time- 
variable labels the subsequent measurements. The aver- 
ages of current powers (noise and third cumulant) have 
been already calculated 0, with respect to the 

quasiprobability (without introducing this notion) and 
measured In experiments, the large measured off- 

set noise plays the role of g in ([3]) preventing paradoxical 



results [12|. We emphasize that in long-time averages 



the quasiprobability becomes a conventional probability 
and reproduces the formula for the usual full counting 
statistics [Isl , also confirmed experimentally [13, [11] . 

We consider a quantum point contact described by 
fermionic operators around the Fermi level [2l|. Each 
operator tpAn with n — (n, a) is denoted by mode num- 
ber n £ {l..A^} and spin orientation a and A = L,R for 
left and right going electrons, respectively. Each mode 
can have its own Fermi velocity Vn and transmission coef- 
ficient Tn (reflection i?„ = 1 — T„). We will assume non- 
interacting electrons and energy- and spin-independent 
transmission through the junction. The Hamiltonian in- 
cluding a voltage bias V reads 

H = ^ dx ^ihvn[ipl^{x)dxipLn{x) - R] 

+qnS{x)[iplf^{x)ipRn{-x) + ^/'ljft(a;)'0Ln(-a;)] 
-e^0(a:)[V;L(a^)^Lft(a:)+^k(^)^i?.ft(^)]} ■ (10) 

The fermionic operators satisfy anticommutation rela- 
tions {'ipa{x),^pb{x')} = and {■ipa{x),tpl{x')} = SabS{x- 
x') for a,b — Ln, Rffi. The transmission coefficients are 
Tn = cofi^{qn/hvn). We apply ^ to the current operator 
/ = J2n^'"n-4'ln{0+)il'Ln{0+) - L ^ R and the density 
matrix p oc ex.-p{—H/kBT). 

To find conditions in which the inequality ([1]) is vi- 
olated, it is enough to consider the case V = Q. Us- 
ing Eq. ([9]) we obtain Sn = hGw{a) and Cap = 
hFGe'^{w{a)+w{P))/2 (Hfii. Here we denote w{oj) = 
u! coth{hu!/2kBT), conductance G = X^n ^^-^"/'"'^ ^^'^ 
Fano factor F = X)n ^"^n/ X^n ^n- I'^ the tunneling 
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FIG. 1: (color online) (a) The minimal value of the con- 
ductance G which satisfies the inequality ((Sjl for the meso- 
scopic junction at zero voltage. The non-classicality occurs 
for G < Gmin • (b) The dependence of inequality parameter 
B [given by ^\ for a the tunnel junction as function of the 
high frequency. Below the line at 1 is the classical regime. Ei- 
ther finite voltage and/or temperature lead to suppression of 
B and result in the necessity to measure at higher frequencies. 



limit, r„ ^ 1, for a finite bias voltage V we only need to 
replace w[uj) with {w{uj + eV/h) + w{uj — eV/h))/2 and 
F = 1. For T = 35mK, 5u = Slo' = 2cj' = 27r • 200MHz, 
a; = 27r • 6GHz, G-^ = 500kl7 and ^ = 0, we get 
B — lA, which contradicts our classical expectation ^ 
and clearly shows that the quasiprobability g must take 
negative values. Generally, the violation occurs for suf- 
ficiently small G < Gmin, as shown in Fig. [T]i. At 
eV ^ and Suj_^ 2uj > Scu' ~ 2a;' > ksT/H, we have 
Gmin = 3y/u;/uj'Fe'^/2h. For larger conductance, one 
can still find a reasonable range of parameters for the vio- 
lation at G~^ = 5kri, as shown in Fig. [TJd. The strongest 
violation occurs at low temperature and voltage but at 
large bandwidth. Unfortunately, the typical experimen- 
tal conductance is with G^^ sa 50Q 29| even larger and 
would require either a; ^ 27r • ITHz or a temperature 
~ ImK. However, we can make the reasonable assump- 
tion that all modes of the junction are independent and 
replace the inequality ([T]) - valid for the whole junction - 
by the same one for a single mode. If we can assume that 
the modes are independent and have similar transmission 
coefficient (T„ ^ 1) we can simply divide G and S by the 
number of modes in ([T]). Note that C enters there linearly 
while S enters quadratically, so effectively we weaken the 
contribution from the second cumulant. For a tunnel 
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junction we can thus replace G by G/Nq where Nq is a 
lower bound of the contributing modes (T„ > 0), which 
must be larger than G-Kh/e^. The results in Fig. [T}) are 
valid for example for = 100 < (/i/2e2)/50f7 « 258. 
In this case, it is necessary to ensure that most of the 
modes contributing to the transport have small trans- 
mission eigenvalues. 

The cumulants are never measured directly. The sec- 
ond cumulant always contains a large offset noise gener- 
ated by detector and amplifier and effectively described 
by A in g. The offset noise can presumably be subtracted 
as it is due to detector's amplifier. The offset noise may 
be smaller in the case of a many-mode tunnel junction 
and one way around is to measure cross correlations by 
different detectors and amplifiers [12|. The fourth cumu- 
lant should also contribute to photon counting statistics 
30l 31 1, but in the limit (5w — > the photon statistics is 
dominated by the second cumulants in ([l}. 

We have shown that the unusual properties of weak 
measurements can be interpreted in terms of a real 
quasiprobability, which can take negative values. Our 
interpretation agrees well with predictions and measure- 
ments of the current fluctuations in mesoscopic junc- 
tions. Its direct confirmation would be the measurement 
of high-frequency fourth-order averages of the current 
though the junction. By a violation of the inequality p]). 
the negativity of the quasiprobability could be directly 
demonstrated. Finally, the separation between detector 
and the system is somewhat arbitrary. One could argue 
that simply adding some noise can restore the positive 
probability. This is why an experimental estimate of the 
detector noise g or the strength of the measurement A is 
also desirable. 
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